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ABSTRACT 

We calculate the projected two point correlation function for samples of luminous and massive galaxies in the 
COMBO- 17 photometric redshift survey, focusing particularly on the amplitude of the correlation function at 
small projected radii and exploring the constraints such measurements can place on the galaxy merger rate. For 
nearly volume-limited samples with 0.4 < z < 0.8, we find that 4± 1 % of luminous Mg < -20 galaxies are in close 
physical pairs (with real space separation of < 30 proper kpc). The corresponding fraction for massive galaxies 
with > 2.5 X 1O'"M0 is 5±1%. Incorporating close pair fractions from the literature, the 2dFGRS and the 
SDSS, we find a fairly rapid evolution of the merger fraction of massive galaxies between z = 0.8 and the present 
day. Assuming that the major merger timescale is of order the dynamical timescale for close massive galaxy 
pairs, we tentatively infer that ~ 50% (70%) of all galaxies with present-day masses M» > 5 x 1O'"M0 (remnants 
of mergers between galaxies with > 2.5 x 10'"Mq) have undergone a major merger since z = 0.8(1): major 
mergers between massive galaxies are a significant driver of galaxy evolution over the last eight billion years. 

Subject headings: galaxies: general — galaxies: interactions — galaxies: evolution 



1. INTRODUCTION 

Galaxy mergers are ubiquitous in a hierarchical universe, 
and are predicted to be an important mode of galaxy growth, 
particularly at early times in cosmic history (e.g., Keres et al. 
2003; Mailer et al. 2006). Mergers may be an important fea- 
ture of the growth of massive early-type galaxies (e.g., Toomre 
& Toomre 1972; Barnes & Hernquist 1996; Khochfar & Burk- 
ert 2003; van Dokkum 2005; Naab, Khochfar, & Burkert 2006; 
Bell et al. 2006). It is possible that mergers may also be an im- 
portant driver of disk galaxy evolution: gas-rich mergers may 
conserve enough angular momentum to form a disk (Robert- 
son et al. 2004). The non-circular gas motions induced by the 
rapidly-changing potential may drive gas inflows, igniting in- 
tense star formation (see, e.g., Sanders & Mirabel 1996, for 
a review), feeding pre-existing supermassive black holes, en- 
hancing AGN activity, and perhaps even driving a galaxy-scale 
superwind, evacuating the galaxy of gas (Springel, di Matteo, 
& Hernquist 2005; di Matteo, Springel, & Hernquist 2005). 

Despite their importance, it has proven challenging to mea- 
sure the rate of major galaxy merging, and its evolution with 
cosmic epoch. The measurement of the galaxy interaction rate 
by counting the incidence of strongly-disturbed galaxies (with 
strong asymmetries, double nuclei, or prominent tidal tails) has 
provided important constraints on merger rate (e.g., Le Fevre 
et al. 2000; Conselice et al. 2003, 2005; Lotz et al. 2006), 
but suffers from uncertainties: minor gas-rich interactions may 
produce much more spectacular results than a major merger 
between two spheroid-dominated galaxies, and the timescales 
over which merger signatures are visible is highly dependent on 
orbits, gas content, and mass ratio. Another powerful method 
for exploring the galaxy merger rate is to measure the inci- 
dence of close pairs of galaxies: it allows access to the prop- 
erties of the progenitors (and therefore, e.g., stellar mass ratio 
of the merger), is straightforward to quantitatively measure, and 
can be modeled using current generations of galaxy formation 



models. Yet, this has proven to be a reasonably challenging 
endeavor: contamination by projection, luminosity boosts by 
interaction-induced star formation, and small number statistics 
are significant challenges and are not easily circumvented. 

In this paper, we present our first attempt at addressing 
this issue using the COMBO- 17 photometric redshift survey 
(§2). This analysis uses a large sample, attempts to correct 
the close pair fraction estimates for projection, and incorpo- 
rates stellar mass estimates, making it highly complementary to 
other important recent attempts at measuring the evolution of 
close pair fraction (e.g., Patton et al. 2002; Lin et al. 2004). 
We split the problem into two aspects. Firstly, we measure 
the fraction of galaxies in r < 30 kpc separation pairs (in real 
space) through analysis of the projected correlation function 
of galaxies from COMBO-17 (§3). This is a well-understood 
and well-posed problem, with a clear and well-constrained out- 
come (§4). We subsequently use estimates of merger timescale 
to explore implications for the merger rate of galaxies: inas- 
much as this part of the analysis makes use of somewhat un- 
certain merger timescales and assumes that all real space close 
pairs will merge, this part of the analysis is much less robust 
(§5). Throughout, we assume ilm = 0.3, r^m + f^A = 1, and 
//o = 70kms-iMpc-i. 

2. THE DATA 

To date, COMBO-17 has surveyed three disjoint - 34' x 33' 
southern and equatorial fields to deep limits in 5 broad and 12 
medium passbands. Using these deep data in conjunction with 
non-evolving galaxy, star, and AGN template spectra, objects 
are classified and redshifts assigned for ^ 99% of the objects 
to a limit of ^ 23.5. Typical galaxy redshift accuracy is 
Sz/(l+z) ^ 0.02 (Wolf et al. 2004), allowing construction of 
^0.1 mag accurate rest-frame colors and absolute magnitudes 
(accounting for distance and A:-correction uncertainties). Astro- 
metric accuracy is ^ 0.1". Owing to reduced depth close to the 
edges of the fields, we discard galaxies < 1' from the image 
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edge. 

As we are concerned with the clustering of galaxies on small 
angular scales, we have tested how the detection of galax- 
ies by COMBO- 17 is affected by having a nearby luminous 
neighbor. We extracted the images of 400 isolated massive 
galaxies (galaxies which were included in our sample) in the 
COMBO- 17 image of the Extended Chandra Deep South, and 
placed these postage stamps of real massive galaxies with a 
distance 0" < r < 10" from another massive galaxy. We then 
put this modified image through the COMBO- 17 object detec- 
tion pipehne, allowing us to determine the fraction of these in- 
serted massive galaxies which were recovered by COMBO-17's 
pipeline as a function of distance from the primary galaxy. We 
found that the detection fraction was independent of distance at 
r > 2"; only at r < 2" (corresponding to ^ 15 kpc at the redshift 
of interest) was there evidence for substantial incompleteness in 
object recovery. Accordingly, we do not use information from 
pairs with separations < 15 kpc in what follows. 

Borch et al. (2006) estimated the stellar mass of galaxies 
in COMBO- 17 using the 17-passband photometry in conjunc- 
tion with a non-evolving template Ubrary derived using the Pe- 
GASE stellar population model (see Fioc & Rocca-Volmerange 
1997, for a description of an earlier version of the model). The 
masses were derived using a Kroupa et al. (1993) stellar IMF; 
the use of a Kroupa (2001) or Chabrier (2003) IMF would have 
yielded similar stellar masses, to within ~ 10%. The redder 
templates have smoothly-varying exponentially-declining star 
formation episodes and a low-level constant star formation rate; 
the bluer templates have a recent burst of star formation super- 
imposed (thus, ongoing tidally-induced bursts of star formation 
are approximately accounted for). Such masses are quantita- 
tively consistent with those derived using a simple color-stellar 
M/L relation (Bell et al. 2003), and comparison of stellar and 
dynamical masses for a few z ~ 1 early-type galaxies yielded 
consistent results to within their combined errors (see Borch 
et al. 2006, for more details). Random stellar mass errors are 
< 0.3 dex on a galaxy-by-galaxy basis, and systematic errors in 
the stellar masses (setting the overall mass scale and its redshift 
evolution) were argued to be at the 0. 1 dex level. 

3. THE METHOD 

As stated in the introduction, our goal is to estimate as accu- 
rately as possible the fraction of galaxies which have a compan- 
ion satisfying our selection criteria within a certain physical dis- 
tance. The accuracy of photometric redshifts is clearly insuffi- 
cient to directly estimate the redshift space correlation function 
of galaxies directly: COMBO-17's typical redshift error trans- 
lates into ~ 200 Mpc along the line-of-sight. Consequently, we 
separate this problem into two parts: construction of a projected 
correlation function, and subsequent de-projection of this cor- 
relation function into a real-space correlation function. 

3.1. Estimating the projected correlation function 

As a first step, we calculate the projected two-point correla- 
tion function of galaxies (Davis & Peebles 1983). The projected 
correlation function w{rp) is the integral of the real space cor- 
relation function ^{r) along the line of sight: 

/oo 
m'^„ + lT^]"^)dlT, (1) 

where is the transverse distance between two galaxies and 
TT is their Une-of-sight separation. It is clear that traditional 



close pair counts are an integral over small projected radii of 
the projected two-point correlation function for galaxies with 
particular properties, multiphed by the space density of such 
galaxies. 

In practice, we estimate wi/p) for various samples of galax- 
ies using the following scheme. We construct a histogram of the 
number of galaxy pairs with given properties (as defined later in 
§4) and | Az| < 0.05 as a function of projected physical separa- 
tion, and for a randomly-distributed mock galaxy sample. The 
mock samples were generated from the real data by bootstrap- 
ping the data multiple times, assigning random positions and 
fields, and applying a small Gaussian redshift offset (ctj = 0.04). 
Masks defining the field edges and the areas around bright stars 
were applied in the same manner to the data and random fields. 
Projected auto-correlations w{rp) were estimated from these 
histograms by constructing the ratio w{rp) = A(DD/RR-1), 
where A is the path length being integrated over, DD is the 
histogram of separations between real galaxies and RR is the 
histogram of separations of mock catalog galaxies. Other for- 
mulations ( rr"*^*^ and DD/DR-1, where DR is the his- 
togram of separations of real and mock galaxies) were verified 
to be equivalent to within the errors. 

One will see that, in order to preserve S/N, we did not inte- 
grate along the entire fine of sight in calculating wirp), rather 
along ±0.05 in redshift path length from the galaxy of inter- 
est. Therefore, we must correct our estimate of w(r^) to ac- 
count for pairs missed because their photometric redshift er- 
rors took them erroneously out of the redshift range being in- 
tegrated over. Extensive comparison with spectroscopic red- 
shifts (Wolf et al. 2004) has shown that a Gaussian with the fol- 
lowing /?-band apparent magnitude-dependent width is an ad- 
equate representation of the photometric redshift errors: ~ 
0.01134[1 + 10"**('"«-2i-5)]i/2. The average a, for the galaxy 
sample of interest was evaluated, and the redshift difference dis- 
tribution of galaxy pairs was taken to be described by a Gaus- 
sian with width (7p.{,r = v^cr,-- Accordingly, the fraction of gen- 
uine galaxy pairs included in the sample is: 

/0.05 1 
e-^'/H^dz, (2) 
0.05 v27r<Tpair 

(i.e. the fraction missed was 1 -/). The estimate of w(rp) was 
then multiphed by 1// to account for the missing pairs, a cor- 
rection of < 30% in all cases. No correction for the integral 
constraint was applied: Phleps et al. (2006) show in their §4.4 
that such corrections are < 0.5 Mpc in w(rp) for 0.4 < z < 0.8 
luminous/massive galaxies in COMBO- 17, which given our fo- 
cus on the strongly correlated smallest scales is a <C 1% correc- 
tion. 

3.2. Estimating the real space correlation function 

Assuming that the real-space and projected correlation func- 
tions can be adequately fit with power laws, the parameters 
of the two fits are intimately related. If ^(r) = (r/ro)~^, then 
w(rp) = Cr^r^'^, where C = v^^To^T^ ' have adopted 

the approach of fitting the w{rp) estimated from the data as a 
function of projected radius, then using the above relations to 
estimate ro and 7. It is also possible to directly estimate ^(r) 
from w(rp) using the Abel integral (Davis & Peebles 1983): we 
choose not to adopt this method by default in this work because 
the direct inversion is rather noise-sensitive, and because we are 
interested in very close pairs r < 30 kpc, and some power-law 
extrapolation of ^(r) would be necessary at any rate to fill in the 
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closest < 15kpc separations (see above). We have confirmed 
that in the best-posed cases where the w(rp) are well-measured 
that i) the correlation functions are well-parameterized by a 
power-law, and ii) inverting w(rp) directly into ^(r) gives very 
similar answers to the power-law parameterization. 

3.3. Estimating the close pair fraction 

Given a good parameterization of the real space correlation 
function of galaxies on small < 30kpc scales, it is straightfor- 
ward to define a real space close pair fraction (as discussed by 
Patton et al. 2000; Masjedi et al. 2006). Recall the definition of 
the real space correlation function: SP = n[l + (,{r)]SV , where 
SP is the probability of a secondary galaxy occupying a volume 
SV a distance r from the galaxy of interest, and n is the space 
density of the secondary galaxies (equal to the space density 
of the primary galaxies for an autocorrelation). Therefore, the 
probability of a galaxy being within a distance rj of another 
galaxy satisfying our selection criteria is: 

Pir<rf)= [ ' n[l+ar)]dV, (3) 

Jo 

«47rn / r^^(r)dr, (4) 
Jo 

given that ^(r) I:^ 1 at all radii of interest for this paper. There- 
fore, parameterizing the real-space correlation function as ^ (r) = 
{r/ro)~^, one obtains: 

. 47rn 

Pir<rf) = - ry.^. (5) 

J — 7 ■' 

It is worth noting that because typically 7 ~ 2, P(r < rf) oc rj 
to first order: i.e., that there is roughly an equal contribution of 
galaxies in each radius bin to the total close pair fraction. 

4. RESULTS 

As stated earlier, our goal is to understand the role of merg- 
ing in driving the evolution of massive galaxies. Therefore, we 
study the close pair fraction of galaxies selected in two ways: 
luminous Mb < -20 galaxies, and massive M* > 2.5 X 1O'"M0 
galaxies. We restrict the sample to galaxies in the redshift in- 
terval 0.4 < z < 0.8. Each sample is ^ 99% complete across 
the entire redshift range of interest, forming a nearly volume- 
Umited sample. The auto-correlation function w(rp) for each 
sample has been derived using the methods outlined above. A 
power law is fit to the w(rp) values for the whole sample, and 
the close pair fraction P(r < r/) derived from this power law fit 
coupled with the measured number density of galaxies meeting 
our selection criteria. Error bars in all quantities were derived 
for each field separately, adopting Monte Carlo errors of 0.1 
dex in luminosity/mass and accounting for counting uncertain- 
ties in the histogram of real galaxy pairs DD. These were then 
combined in quadrature and divided by VC-^fieid- 1) = V2. It 
is important to note for this paper we have used proper coordi- 
nates to calculate the correlation functions and space densities. 

The results are shown in Figs. 1 and 2. One can see that the 
power-law parameterization to w(rp) is an acceptable descrip- 
tion of the data on sub-Mpc scales, and that w(rp) is measured 
with interesting accuracy even to ~ 20 kpc scales (allowing rea- 
sonably robust measurement of P(r < rf)); we showed earUer 
that object detection failed for separations < 15 kpc, therefore 
power law fits to w(rp) are determined only for 15 < r/kpc < 
1000, and are extrapolated inwards. The resulting close pair 
fractions Pir < rf), adopting rf = 30 proper kpc, are shown in 
Fig. 2. 
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Fig. 1. — Grey: Projected correlation function »(?■,,) for tlie 0.4 < z< 0.8 
volume-limited sample of luminous Mb < -20 galaxies. The power law fit 
to w(rp) is overplotted, and parameters given in Table 1. The vertical dashed 
line at 15 kpc shows the radius within which COMBO-17's object detection 
pipeline no longer reliably separates nearly equal-luminosity close galaxy 
pairs: this corresponds to the radius at which the correlation function starts to 
deviate strongly from a power law. Black: The projected correlation function 
of the volume-limited 0.4 < z < 0.8 sample of massive M* > 2.5 x IO^^Mq 
galaxies. 



Table 1 

Real space correlation function parameters and close 

PAIR fractions 



Sample ro/Mpc 7 nMpc^ P(r < rf) 

Auto-correlations 

(T) 3.1±0.5 1.93±0.08 00162 0.04±0.01 

(2) 3.6±0.5 2.02 ±0.07 0.0091 0.05±0.01 

Cross-correlations 

(3) 3.2 ±0.5 1.89 ±0.09 0.0136 0.028 ±0.005 

(4) 3.7 ±0.5 1.98 ±0.08 0.0061 0.029 ±0.005 



Sample (1) is a luminosity and volume-Umited sample with 
Mb < -20 and 0.4 < z < 0.8 

Sample (2) is a stellar mass and volume-limited sample with 

> 2.5 X IO^^Mq and 0.4 < z < 0.8 
Sample (3) has primary galaxies with Mb < -20.3 and 
0.4 < z < 0.8; secondary galaxies must be 1.2 mag fainter or 
less than their primary galaxy 

Sample (4) has primary galaxies with M* > 3 x lO^^M©; 
secondary galaxies must be 1/3 or more of their primary 
galaxy's mass to be counted as a pair 

For samples (3) and (4), the space density of secondary 
galaxies is used to determine the pair fraction 
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Many previous studies have attempted to estimate the ma- 
jor merger rate, i.e., the rate of galaxies merging with galaxies 
with mass ratios between 1:1 and 3:1. Accordingly, we have 
also estimated the cross-correlations between galaxies in our 
sample and potential major merger partners: for the luminous 
galaxy samples, between Mb < -20.3 galaxies and galaxies be- 
tween and 1.2 mag fainter than the primary galaxy; and for 
the massive galaxy samples between > 3 x 10^^ Mq galax- 
ies and galaxies with between 1/3 and the same mass as the 
primary. The samples are close to volume limited in each case; 
the faintest red sequence galaxies are missing at z > 0.6 in the 
secondaries (the blue cloud secondary galaxies are complete 
at all redshifts), leading to a < 7% incompleteness in the sec- 
ondary sample; we do not correct for this incompleteness. The 
number density of secondary galaxies is not trivial to calculate 
as each primary galaxy has a different set of secondary galax- 
ies. We estimated the 'characteristic' number density of sec- 
ondary galaxies by evaluating the number density of potential 
secondary galaxies for each primary galaxy, and then averag- 
ing these number densities. The cross-correlation parameters 
are very similar to those for the autocorrelations, with the mod- 
est differences between the values of P(r < rj) for the auto- 
and cross-correlations driven by the differences in the number 
densities. We give parameters of power-law fits to the w{rp) 
estimates in Table 1. 

5. DISCUSSION 
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Fig. 2. — The close pair fraction (with real space separations < 30kpc) 
of luminous {Mb < -20; left panel) and massive (M» > 2.5 X IO^^Mq; right 
panel) galaxies. The z ~ 0.6 data points are from the present work, and the 
z ~ 0.1 data points are estimated (roughly) from the 2dFGRS (left) and SDSS 
(right). In the case of the 2dFGRS, two estimates are shown: the lower estimate 
is for Mb < -20 galaxies, while the upper estimate is for galaxies with Mb < 
-19 A — the plausible descendants of Mb < -20 galaxies at z = 0.6. 



5.1. Comparison with local pair fractions derived using the 

same method 

Using published studies, one can use the methodology pre- 
sented in this paper to estimate the local pair fraction. For 
luminous galaxies, we use results from the Two degree Field 
Galaxy Redshift Survey (2dFGRS; Colless et al. 2001). Nor- 
berg et al. (2002a) presented the two-point correlation func- 
tion for galaxies with Mi,^ -51ogjQ/iioo < -19.5 (correspond- 
ing to Mb < -20.1 adopting bj B-0.15 from Norberg et 



al. and converting to /iioo = 0.7). The values of vq and 7 are 
3 .5 ± 0.4 Mpc and 1 .8 ± 0. 1 respectively, and were defined only 
outside IOO/iJqq kpc; the inwards extrapolation required for this 
analysis is unconstrained and therefore this estimate should 
be regarded with due caution. Adopting the luminosity func- 
tion from Norberg et al. (2002a), a density of galaxies with 
Mb < -20 of 0.0019(2)Mpc~^ was derived, giving a pair frac- 
tion estimate of 0.003(1). An unavoidable compUcation is that 
of the fading of galaxies towards lower redshift as their stellar 
populations fade; thus, a sample of somewhat fainter galaxies 
at lower redshift may be a better conceptual match to the pop- 
ulation of distant Mb < -20 galaxies. Adopting a somewhat 
fainter cut of Mb{z-q) < -19.4 (corresponding to 1 mag fading 
per unit redshift in the rest-frame B-band), one derives instead 
a largely unchanged clustering signal (Norberg et al. 2002a), a 
larger density of 0.0039(4) Mpc"^, and a larger pair fraction of 
0.006(2). Two points become clear: two-point correlation func- 
tion parameters are substantially less sensitive to limiting depth 
than the number density or close pair fraction; and, the fading 
of stellar populations is a considerable complicating factor in 
studying the evolution of close pair statistics for luminosity- 
selected samples. 

Stellar mass-limited samples largely overcome the last of 
these two challenges. A large number of galaxies in the main 
galaxy sample of the Sloan Digital Sky Survey (SDSS; York 
et al. 2000) have estimates of their stellar mass. We use the 
two-point correlation function of mass-limited samples (Li et 
al. 2006) and the space density of massive galaxies (Bell et 
al. 2003). Rough fits to Li et al.'s 10.2 < logM^^/M© < 10.7 
and 10.7 < log M^/M© < 11.2 samples gave correlation func- 
tion parameters ro ^ 4.2 ±0.4 and 7 ^ 1.85 ±0.05: given a 
number density of > 2.5 x lO'^M© galaxies of 0.0040(4) 
galaxies/Mpc^ from the SDSS stellar mass function of Bell et 
al. (2003), we estimate a P(r < 30kpc) ~ 0.011(5). Again, the 
correlation function was defined only outside lOO/zJookpc; for 
both the SDSS and the 2dFGRS a direct redetermination of the 
pair fractions from the data would be preferred to these rough 
estimates. The average redshift of the galaxy samples in both 
cases is (z) 0.1. 

These estimates are included in Fig. 2 as the z ~ 0.1 esti- 
mates. It is worth noting in both cases that there is clear evi- 
dence for a dramatically-reduced fraction of galaxies in close 
(< 30 kpc separation) physical pairs at the present day, com- 
pared to the pair fraction at z ~ 0.6. This will be discussed in 
more detail later. 

5.2. The relationship between true close pairs and projected 

close pairs 

In order to compare our measurements to others in the lit- 
erature, it is necessary to explore the relationship between the 
fraction of galaxies in real space close pairs and the projected 
close pair galaxy fraction. The projected space close pair frac- 
tion is the real space close pair fraction plus a contribution from 
more distant galaxies along the line of sight^ Specifically, the 
projected close pair fraction: 

P\rp<rf) = P{r<rf)+ I A'Kr^n[\+^{,r)\\_2/T:sm\rf /r)fdr, 

(6) 

following previous notation. In our particular case, r/ = 30 kpc, 
and using the observed luminous/massive galaxy correlation 

'it is worth remembering that the projection of very distant fore- and back- 
ground galaxies has been automatically removed when calculating wirp). 
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functions and number densities as input, we find that 65% 
of projected close luminous pairs have real space separations 
< 30kpc; the corresponding fraction for massive galaxies is 
69%. It is important to note that this 'contamination' is with 
galaxies which are correlated with the host (i.e., those which are 
primarily nearby with 30 < r/kpc < 1000), and would likely 
have very similar redshifts to the primary galaxy (i.e., much 
of contamination is suffered by spectroscopic close pair sam- 
ples). It is important to note that the exact fraction depends on 
the detailed form of the correlation function and should not be 
blindly adopted by workers using rather different sample cuts: 
in particular, our estimate is sUghtly higher than the estimate 
of ^ 50% from Patton et al. (2000), which was derived using 
a very similar approach with a less clustered lower-luminosity 
parent sample. 

5.3. Comparison with published merger fraction 
determinations 
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Fig. 3. — The close pair fraction (with real space separations < 30kpc) 
of luminous {Mb < -20; left panel) and massive (M, > 2.5 x 10"*Mq ; right 
panel) galaxies. Left: Gray data points show the merger fraction of galaxies 
with Mb < -20, taken from a variety of soiffces (see §5.3 for details). The 
Somerville et al. model prediction of the major merger fraction for Mb < -20 
galaxies is shown as a solid gray line for the redshift range 0.1 < z < 1.1. 
The z ~ 0.1 data point is derived from the 2dFGRS. Right: Gray solid circles 
show the merger fraction of galaxies with M* > IO^'^Mq from Conselice et al. 
(2003), and the Somerville et al. model prediction of the major merger rate of 
massive galaxies is shown as a solid gray line. The merger fraction predicted 
by Mailer et al. (2006) is shown by the grey dotted line. The z ~ 0.1 black data 
point is derived from the SDSS, and the gray diamond is taken from Xu, Sun 
& He (2004). 

In Fig. 3, we have attempted to show close pair or merger 
fractions from a variety of published works, to compare with 
our determinations. Such an exercise is not trivial, as very dif- 
ferent methods were used in many cases. 

5.3.1. The Mb < -20 sample 

The vast majority of previous determinations of close pair 
fractions have been derived for luminous Mb < -20 galax- 
ies. At z < 0.3, we show results for Patton et al. (2002) 
and de Propris et al. (2005) as open circles. Patton et al. 
(2002) adopts a radius range 5/zjqq < r < 20/z7QQkpc; adopt- 
ing Hq = 70kms"' Mpc~' this becomes 7 < r < 30kpc (i.e., 
some 1/4 of the galaxy pairs are not counted). Patton et al.'s 



default magnitude range is -21 < MB-51ogiQAioo < -18; cor- 
responding to -22 <Mb< -19 for Hq = 70kms"> Mpc"^ Us- 
ing their table 3, one can convert their results to a narrower 
range in absolute magnitude -21 < MB-51ogjQ/!ioo < -19, or 
-22 < Mb < -20 in our units, by dividing by 2.5. Accordingly, 
we adjust Patton et al.'s values upwards by 4/3 (to account for 
r < 7kpc pairs) and downwards by a factor of 2.5 (to account 
for the magnitude range), de Propris et al. (2005) adopt a mag- 
nitude range of -22 < -5 logjQ hioo < -19, corresponding to 
-23 < Mb ^ -20 for our choice of Hq, thus the only correction 
appUed is 4/3, to account for missed r < 7kpc pairs. 

At z > 0.3, we show the fraction of Mb < -20 galaxies in 
close pairs with < 30kpc separation (their < 20h'^l^kpc val- 
ues) taken from Le Fevre et al. (2000, open diamonds); we have 
adjusted the values downwards to 65% of their original values 
to account for projection within galaxy groups (as discussed 
above, as opposed to projection of random galaxies along the 
line of sight, which Le Fevre et al. corrected for already). We 
do not show values from Bundy et al. (2004), whose optical 
pair statistics agree well with Le Fevre et al.'s values but who 
argue (based on near-infrared data) that many of the appar- 
ently luminous pairs are in fact minor mergers which have been 
boosted in rest-frame B-band luminosity by enhanced star for- 
mation^. We attempt also to include the estimates of Lin et al. 
(2004). They adopt Hq = 70kms"> Mpc"> for the purposes of 
quoting A:-corrected magnitudes, thus their evolution-corrected 
— 21 < Mg < -19 sample corresponds roughly to -22 < Mb ^ 
—20, remembering that the evolution correction is roughly 1 
magnitude per unit redshift. They quote their pair fractions 
in terms of hioo = 1 (L- Lin, 2006, priv. comm.) thus their 
10 < r//!"'kpc < 30 bin corresponds to 15 < r < 42kpc adopt- 
ing our H(). Since we find that P{r < rj) oc rj, their 27kpc of 
coverage for their pair fraction should be approximately equal 
to the P(r < rj) which we would calculate within 30kpc. We 
do, however, apply a correction of 0.65 to their measurements, 
to account for projection at small radii (following §5.2). It is 
clear that the COMBO- 17 0.4 < z < 0.8 estimate is quantita- 
tively consistent with these estimates, to within the combined 
uncertainties, with the advantage of robust projection correc- 
tion, a volume-limited galaxy sample, large sample size, and 
therefore highly competitive errors. 

Owing to the difficulty in extracting the properties of pro- 
genitors from a morphologically-classified ongoing merger, we 
elect not to compare explicitly with morphologically-derived 
merger fractions from Le Fevre et al. (2000), Conselice et al. 
(2003), Cassata et al. (2005), and Lotz et al. (2006). Lotz et 
al. (2006) compare their merger fractions to the other morpho- 
logical studies, finding consistent results to within their com- 
bined error bars. They further compare their results with Lin et 
al.'s and Patton et al.'s results for pair fraction evolution, finding 
overall consistency in both the inferred zero point and redshift 
evolution. 

5.3.2. The > 2.5 x IO'^Mq sample 

To the best of our knowledge, there is only one pubUshed 
estimate of the massive galaxy merger fraction at intermedi- 
ate redshift, for > lO'^M© galaxies from Conselice et al. 
(2003). This estimate is of limited applicability: not only is 

^We address this soitfce of concern through the analysis of the stellar mass- 
limited sample, but find a roughly equal pair fraction. On one hand, many pairs 
of limiinous galaxies are indeed minor mergers by mass; yet, on the other hand, 
there are a large number of red and lower-luminosity galaxies which are missed 
by the luminous galaxy criterion which make it into a mass-limited sample. 
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the mass limit substantially lower, but it is derived from on- 
going mergers, making it almost impossible to fairly compare 
with our pair-based estimate. Furthermore, the small sample 
size weakens their constraints on the merger fraction; in partic- 
ular, they can only place an upper limit at intermediate redshift. 
Nonetheless, their constraints are consistent with ours, to within 
the combined uncertainties. 

Xu, Sun & He (2004) presented an analysis of ^T-selected 
galaxy pairs taken from a combined 2MASS/2dFGRS sample 
(Cole et al. 2001). Converting their results to our value of Hq 
and stellar IMF, extrapolating their 7.5 < r < 30kpc results to 
r < 30kpc, and accounting for the ~ 30% contamination of 
their pair sample with 'group interlopers', we find a pair frac- 
tion of 1%±0.5% for galaxies with > 2.5 x 10'°Mq, in 
excellent agreement with the SDSS determination. 

5.4. The relationship between close pair fraction and merger 

fraction 

In order to compare with galaxy formation models and to 
explore the implications of these and other close pair fraction 
determinations for galaxy merger rate, it is important to discuss 
the relationship between galaxy close pair fraction and merger 
fraction (see, e.g., Lotz et al. 2006, for a recent discussion). 

Let us, for the sake of argument, take the case of galaxy 
mass limited samples. In this work, we derive the close pair 

fraction /pair = A/galpair/A'M.>2.5xl()i»Mo, whcrC A^gaLpair IS the 

number of galaxies in pairs with M, > 2.5 x IO^^Mq, and 
^M,>2.5xioi»Mo is the number of massive galaxies, per given 
volume. If one wanted to relate this close pair fraction to a 
merger fraction, i.e., the fraction of galaxies which have been 
created by mergers of galaxies in the pair sample, one notes 
that /) two galaxies in pairs merge into one merger remnant, 
ii) the merger remnants are higher mass, and the timescales 
of being a recognizable merger remnant may differ from the 
close pair timescale. In this example above, the merger frac- 
tion most directly related to the above close pair fraction is: 
/merg = A^merg/A^M. >5 X ioi»Mq ; i-^-, the uumbcr of ucwly-crcated 
merger remnants with M* > 5 x XO^^Mq (the remnants of merg- 
ers between galaxies with > 2.5 x lO'^M©) divided by the 
number density of M* > 5 x lO'^M© galaxies. The number of 
pairs should be then be related to the number of newly-created 

remnants. A^gal.pair ~ 2A^merg^pair/'^merg) 

where Tpair/rmerg is the 
ratio of the timescales over which a pair enters one's close pair 
sample vs. the timescale over which a merger renmant is rec- 
ognizably disturbed, and the factor of two accounts for the fact 
that a galaxy pair merges to form a single remnant. 
Let us make this example more concrete. 

/pair = A'gal,pair/A^M,>2.5xlO"'M0 (7) 

2A^merg'n^air/'^merg ,o\ 

~ — (8) 

■'Vm,>2.5x1O1»M0 
_ 2A'niergTpair/''merg 

- — 77 ) \") 

Q;^'M.>5x1O1°M0 

where a is the ratio in number density between galaxies with 
> 2.5 X 10'"Mq and > 5 X 10'°Mq. This ratio a ~ 2 
in this case (as directly measured from the dataset), giving 
/pair ^ /merg ''pail- /imerg- Thus, the fact that two galaxics in pairs 
merge to form only a single remnant (entering the numerator of 
the fractions) is canceled out by the factor of two different num- 
ber densities (in our particular case) between the pair parent 
population and the plausible newly-created remnant population 



(entering in the denominator of the merger fractions). Many 
previous analyses neglect or underestimate this difference in the 
number density of the parent population from which the pairs 
are drawn vs. the number density of the higher mass merger 
remnants; this is one of the main contributors to our higher in- 
ferred merger rate than those estimated by Lin et al. (2004, see 
§5.6 for further discussion). A similar argument, with similar 
outcome, applies to the luminous galaxy sample. We will use 
this argument in what follows. 

5.5. Comparison with galaxy formation models 

In this section, we compare model major merger remnant 
fractions with our measurements of close pair fraction (see 
Berrier et al. 2006, for a detailed discussion of model insights 
into the meaning and evolution of close pair fraction). In what 
follows, we adopt r^erg = Tpair = 0.4 Gyr, following the es- 
timated pair timescale calculated in the next section. Thus, 
/merg ~ /pair, as the timcscalcs are defined to be equal, as long 
as we choose to explore the fraction of galaxies recently created 
in galaxy mergers with > 5 X 10'OMq and/or Mb < -20.75. 

In both panels of Fig. 3, we show major merger fractions 
taken from an updated version of the Somerville et al. semi- 
analytic galaxy formation model (see Somerville & Primack 
1999; Somerville, Primack, & Faber 2001, for a description of 
the basic model ingredients). The model includes standard pre- 
scriptions for gas cooling, feedback, and dust extinction and 
star formation. Quiescent star formation is parameterized as in 
de Lucia et al. (2004), and bursts of star formation are triggered 
by major and minor mergers, based on results from hydrody- 
namic simulations of merging galaxies (Cox 2004). Feedback 
from AGN is not included in this model, although we found that 
merger rates calculated from a model including AGN feedback 
were similar to those presented here. The model reproduces 
reasonably well the evolution of the luminosity and stellar mass 
function of galaxies in the interval < z < 1. This particular 
model uses Monte Carlo reaUzations of dark matter merger his- 
tories based on the analytic Extended Press-Schechter formal- 
ism, supplemented with standard prescriptions for dynamical 
friction, and therefore we lack detailed information about the 
spatial location of galaxies within their dark matter halos. As 
a result, we cannot directly compute close pair fractions from 
these simulations; rather, we compute the fraction of galaxies 
which have undergone a recent merger and compare this with 
the observational estimates. 

In the right-hand panel of Fig. 3, we show also the close 
pair fraction inferred from the merger rates of massive galaxies 
from a SPH galaxy formation model (Mailer et al. 2006, dotted 
line), under the same timescale and merger vs. pair assump- 
tions as adopted above. Mailer et al. (2006) present the fraction 
of mergers per Gyr for galaxies with M* > 6.4 x lO'^M© with 
mass ratios less than 2:1 in their Fig. 3: they argue that this 
is equivalent to a mass cut of ^ 2.5 x IO'^'Mq, as their model 
dramatically over-produces stellar mass by a factor of 2.75, and 
as a stop-gap measure they suggest division of the mass cut by 
that factor when comparing with data. This mass limit is differ- 
ent from the merger remnant mass limit of > 5 x IO^^Mq; 
accordingly, we treat this comparison as more qualitative than 
quantitative. The merger rates are increased by a factor of 1.7 
to transform from 2:1 to a 3:1 threshold (following their Fig. 
6). These merger fractions are somewhat higher than observed, 
and show a steeper redshift dependence than the Somerville et 
al. model estimates. 
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Recently, Berrier et al. (2006) discussed the evolution of the 
close pair fraction of luminous galaxies, finding overall con- 
sistency at the factor of two level with the Lin et al. (2004) 
measurements; given that our measurements are consistent with 
the measured pair fractions from Lin et al. (2004), we would 
expect that the model close pair fraction from Berrier et al. 
(2006) would be reasonably consistent with both our data and 
the Somerville et al. model. 

While important discrepancies remain between different mod- 
els, different datasets, and between the data and models, it is 
nonetheless obvious that there is an overall qualitative consis- 
tency between our best attempts at observationally constraining 
the merger/close pair fraction of luminous and massive galaxies 
and our present understanding of galaxy formation and assem- 
bly in a ACDM universe. Given that the ongoing assembly of 
massive galaxies is a key (and unavoidable) feature of galaxy 
evolution in such a cosmology, it is encouraging that there is 
a decent qualitative agreement between the models and data at 
this stage. 

5.6. Musings on the merger rate of galaxies 

In order to convert close pair fractions into merger rates, a 
timescale over which a close pair of nearly-equal mass galaxies 
will merge is required. The estimation of such timescales far 
from straightforward and is the topic of much ongoing work: 
the mix of orbital parameters will lead to a distribution of 
timescales, and the effects of e.g., dynamical friction and fly- 
bys in dense environments are poorly-understood. Here, for 
illustrative purposes, we take a highly simplistic approach and 
assume that the merger of two nearly-equal mass galaxies takes 
roughly one orbital timescale forb ~ 2iTr/\Aa ^ Ar/a, noting 
that the circular velocity of a galaxy is ~ 1 Aa, where a is the 
velocity dispersion of the galaxy in question. This estimate is 
rather similar to those presented by, e.g., Patton et al. (2002) or 
Lin et al. (2004), and compares favorably to timescales derived 
from the Naab, Khochfar, & Burkert (2006) dry galaxy merger 
simulations. We adopt a typical velocity dispersion of a lumi- 
nous/massive galaxy of ~ 150kms~\ and we adopt the typical 
radius of a galaxy in the rj < 30kpc sample of r ^ 15 kpc. On 
this basis, we estimate a merger timescale of ~ 0.4 Gyr; uncer- 
tainties in this timescale are at least a factor of two^. RecalUng 
that /merg ^ /pair whcu cqual timescales are chosen (from §5.4) 
and when the pair fraction from lower luminosity /mass samples 
is related to the merger fraction of the more massive remnant 
galaxies, the merger rate per Gyr is P(r < 30kpc)/0.4. 

Given this rough estimate for merger rate, one can derive the 
rate of creation of galaxies with > 5 x lO'^Mo through ma- 
jor mergers using the estimated close pair fractions of galax- 
ies with > 2.5 X 1O'°M0 at z - 0.1 from SDSS coupled 
with the z 0.6 determination from COMBO- 17. A power- 
law fit to these data points yields: log/ ~ -2.1 ± 0.2-1- [3.9 ± 
0.51og(l +z)], where / is the close pair fraction. Using that 
the creation rate of merger remnants with > 5 x 10 "^M© 
is //0.4, and integrating as a function of cosmic time, one 
derives an average of ~ 0.5 (0.7) galaxy mergers per present- 
day > 5 X 10 '"M© galaxy since z = 0.8(1). The corre- 
al! is worth noting that if dynamical friction timescale arguments were used, 
one would derive timescales for r ~ 15 kpc of ~ 0.3 Gyr, following the discus- 
sion of dynamical friction in Binney & Tremaine (1987). While it is true that 
the Chandresekar formulation of dynamical friction should not strictly apply 
in the case of a merger between two massive galaxies, it is nonetheless en- 
couraging that the orbital timescale and dynamical timescale arguments yield 
approximately equal timescales. 



sponding numbers for galaxies with Mb < -20.75, using all 
available observational determinations and weighting by the 

uncertainties, are: log/ 1.9 ± 0.2 + |2.6 ± 0.81og(l -Hz)], 

and ~ 0.5(0.6) mergers since z = 0.8(1). This determination 
is very similar to the recent determination of Lotz et al. (2006, 
0.3-0.7 mergers per luminous galaxy)'*. While it is clear that 
improvements in both the observational determination of close 
pair fraction evolution with redshift and work towards robust 
estimates of merger timescale are of critical importance, this 
analysis strongly suggests that galaxy mergers between lumi- 
nous/massive galaxies are a significant feature of the last eight 
bilUon years of galaxy evolution. 
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